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Abstract 

As a part of the electrothermal propulsion 
plume research program at the NASA Lewis Research 
Center, efforts have been initiated to analyti- 
cally and experimentally investigate the plumes of 
resistojet thrusters. The method of G.A. Simons 
for the prediction of rocket exhaust plumes is 
developed for the resistojet. Modifications are 
made to the source flow equations to account for 
the increased effects of the relatively large 
nozzle boundary layer. Additionally, preliminary 
mass flux measurements of a laboratory resistojet 
using CO? propellant at 298 K have been obtained 
with a cryogenically cooled quartz crystal micro- 
balance (QCM). There is qualitative agreement 
between analysis and experiment, at least in terms 
of the overall number density shape functions in 
the forward flux region. 

This paper also presents estimates of the 
level of environmental contamination produced by a 
Space Station resistojet during typical propulsion 
operations. Predicted number density profiles 
were used to calculate molecular number column 
density (NCD) as a function of line-of-sight angle 
for an observer located at the astronomical obser- 
vation site of the Space Station. For the case of 
a resistojet using H 2 O propellant, calculated 
NCD values rise above the proposed requirements 
for infrared (IR) species in a restricted viewing 
region representing 20 percent of the total 4ir 
steradian space. In the case of a resistojet 
operating on a benign mixture (e.g., noninfrared 
active) of N2» H 2 , He and Ar, quiescent limits on 
NCD are not exceeded until the observer views 
along line-of-sight angles which pass near, or 
look into, the resistojet source. 

List of Symbol s 

A Simons' plume normalization constant 

A e nozzle exit area, m 2 

A* throat area, m 2 

a 0 sonic velocity at T 0 , m/sec 

C capture coefficient 

f ( e ) Simons* functional variation of density 
M Mach number 

m mass flow rate, kg/sec 

MW molecular weight 

P gas pressure. Pa 


R radial distance from centerline or nozzle 
radius, m 

R universal gas constant 

r spherical radius, m 

S spherical surface area in space, m 2 

T gas temperature, K 

U gas velocity, m/sec 

[i t limiting velocity of the gas, m/sec 

average limiting velocity of the gas in the 
boundary layer, m/sec 

a ratio between U t and U A 

B plume parameter 

y specific heat ratio 

5 boundary layer thickness, m 

a* displacement thickness, m 

0 angle from plume centerline, radians or 
degrees 

e 0 angle corresponding to the streamline at the 
edge of the boundary layer, radians or 
degrees 

ei limiting turning angle of the gas at the 
nozzle exit, radians or degrees 

e„ the value of ei for inviscid, supersonic 
flow, radians or degrees 

y gas viscosity (as a function of temperature) , 
N sec/m 2 

v Prandtl-Meyer (P-M) turning angle based on 
M e , degrees 

p gas density, g/cm 2 or number /cm 2 
4 > nozzle half angle, degrees 
Sub- or Superscript Notation 
* denotes sonic conditions 

e denotes conditions at the nozzle exit 

n denotes limiting 

max denotes maximum 
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o denotes stagnation conditions 

red denotes reduced effecti ve area 

Introduction 

The resistojet thruster has been baselined 
for Space Station auxiliary propulsion applica- 
tions (Level B Space Station Control Board Direc- 
tive). It offers low thrust propulsion for drag 
make-up and also the capability for propulsive, or 
nonpropul si ve, venting of waste gases. Efforts are 
currently directed toward demonstrating the tech- 
nology readiness of a long life, multipropellant 
resistojet for the Space Station. 1-3 

The resistojet concept is relatively simple, 
as illustrated in Fig. 1. A propellant is passed 
through a heat exchanger and then expanded through 
a nozzle to produce thrust. The actual mechanism 
of heat exchange may be one of convection, conduc- 
tion, and/or radiation. Potential resistojet pro- 
pellants include inert gases, water vapor, carbon 
dioxide, oxygen, cabin air, methane, and hydrogen. 
Figure 2 presents a photo of an engineering model 
resistojet designed and fabricated by a Rocketdyne/ 
Technion industrial team under contract to NASA 
Lewis. The design goals are a 10 000 hr lifetime 
and a multipropellant capability. Reference 3 
provides further details of the design and fabri- 
cation of this thruster. 

Flow field definition of Space Station pro- 
pulsion systems and other sources of effluents is 
of fundamental importance. Of primary concern 
are the impacts of contamination, both in terms 
of mass deposition and also environmental contam- 
ination. Consequently, designers of the Station 
core and potential users require definition of the 
resistojet exhaust in order to assess the plume 
impacts to Space Station science and technology 
activities. Other impacts, such as thrust losses 
or thermal loading caused by plume impingement, 
may also be of concern depending upon Space 
Station architecture. 

In order to predict possible contamination 
impacts, the induced environment produced by a 
resistojet propulsion module needs to be modeled 
analytically and, where possible, verified exper- 
imentally. Both tasks are difficult ones. A 
detailed numerical treatment of the plume expan- 
sion begins prior to the exit of the resistojet 
nozzle. Although computer codes to analyze the 
nozzle flow of higher thrust (222 N) propulsion 
devices do exist, the highly viscous nature of 
resistojets, combined with small nozzle geometries, 
precludes the immediate adaptation of these codes 
for the present analysis. Consequently, in order 
to achieve reasonable analytical results in as 
timely a manner as possible, it is instructive to 
begin with source flow analyses of the resistojet 
exhaust. 

G.A. Simons 4 has presented a widely used 
technique for the prediction of a rocket plume 
structure. His method is an extension of earlier 
work by Boynton^ who used finite difference 
numerical techniques to analyze an inviscid plume 
expansion. Viscous boundary layer effects were 
accounted for in the initial conditions. This 
paper presents an analysis of the exhaust flow 
field of a resistojet using Simons' method with 
slight modifications. The method accounts for the 


expansion of the supersonic boundary layer of the 
nozzle, along with the inviscid core. Major plume 
properties are expressed in terms of conditions at 
the nozzle exit and in terms of the boundary layer 
thickness. 

In a parallel effort with the analysis, exper- 
imental mass flux data have been obtained on a 
laboratory resistojet using unheated CC>2 as the 
propellant. A partial mapping of both the forward 
and back flux regions of the exhaust was conducted 
using a temperature-compensated, cryogenically 
cooled quartz crystal microbalance. Preliminary 
mass flux data are then compared with plume den- 
sity calculations. Details of the measurement 
technique and experimental hardware are described. 

Lastly, this paper presents estimates of the 
level of environmental contamination caused by a 
Space Station resistojet during typical propulsion 
operations. The assessment is made for a resisto- 
jet using H 2 O propellant and also a benign gas 
mixture as the propellant. Predicted levels are 
compared with proposed Space Station quiescent 
requirements. 

Analyses 

Modification of Simons' Method 

This section summarizes the development of 
Simons' technique for plume flow field calcula- 
tions 4 and the appropriate modifications appli- 
cable to resistojet exhaust plumes. A detailed 
discussion of the rationale and methodology of 
the following analysis may be found in Ref. 6. 

In the far field, rocket nozzle flow may be 
modeled as a point source. Using continuity, the 
mass flow in the plume may be written in terms of 
conditions at the nozzle throat, as 

p*u* UR* 2 ) = JJp(r,6)\} dS (1) 

It is assumed that the local density p(r,e) can be 
be related to its centerline value p(r,e = 0) by a 
function f(e) such that 

p(r,e) = p(r,8 = 0) f(e) (2) 

Substitution into Eq. (1) gives, upon integration. 


pV (irR* 2 ) = 2wr 2 p(r,e = 0)U t I f(e)sin e de 

•'o 

(3) 

where is the limiting turning angle of the 

gas. The gas velocity U is removed from the 
integrand and becomes the limiting gas velocity 
which is independent of r and e. 

Simons then defines the plume normalization 
constant A, as 
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where = [ Y + 1/ Y - l] U*, based on an 
isentropic expansion to Mach infinity. 

Hence, the local plume density is given by 



The plume structure based on this far field 
approximation is illustrated in Fig. 3. As shown, 
the plume is divided into two, axi symmetric 
regions: the region which originates from the 

inviscid core of the nozzle, and the region which 
originates from the supersonic boundary layer. 

The angle e 0 corresponds to the streamline at 
the edge of the boundary layer. The boundary 
layer is treated as supersonic in this analysis. 
(It should be noted that preliminary results with 
a three-dimensional, Navier-Stokes numerical solu- 
tion of low Reynolds number nozzle flow have 
indicated that a portion of the nozzle may contain 
subsonic flow. However, the code is not suffi- 
ciently developed for incorporation into the pres- 
ent analysis which is meant to serve as a starting 
point for analyses of low thrust, electrothermal 
thrusters. ) 


= 2ttU / pr 2 sin e de (9) 
6 o 

where 5 is the boundary layer thickness. It is 
assumed that the displacement thickness is small 
in comoarison to 5 . The above equation differs 
from Simons' basic equation only in the retention 
of the a 2 term. Simons neglected this higher 
ordered term because he was investigating rocket 
nozzle flow with very small boundary layers 
(i.e., e 0 -*» ej. In resistojet plume analysis 
this term is not negligible because, in many 
cases, the boundary layer may make up a consider- 
able portion of the nozzle volume. As a result of 
the velocity profile in the nozzle boundary layer, 
it should be noted that, in the plume boundary 
layer region, there is more than_one limiting 
velocity such that U becomes^ U* upon its 
removal from the integral. (U* * a where 
0.5 < a < 1.0 for the supersonic boundary layer.) 

Initially, the boundary layer mass flow may 
be treated as inviscid in order to find the value 
for the streamline denoted by the angle e 0 . 

Thus, Eq. (6) is used in Eq. (9) to give 


Simons determined that the appropriate selec- 
tions of f(e) are: 

2/y-l 

for 0 < 6 < 6 (6) 

— — o 

and 


f(e) 





Again, an approximation for sin e is used to 
simplify evaluation of the integral. This gives 
the following expression for e 0 . 



f ( e ) = f (e Q ) exp[- & (e - e Q )] for e Q _< e < (7) 

where e* is the value of for an inviscid, 
supersonic flow. In order to define the density 
contours as given by Eqs. (5) to (7), the values 
of the plume constant A, the parameter 3 , and 
the angle e 0 must first be calculated. It will 
be shown that 6g and 8 are functions of the 
nozzle exit conditions and boundary layer thickness. 

The plume normalization constant requires 
further consideration. If the nozzle mass flow 
were to undergo an inviscid expansion, the dis- 
tribution in the plume would follow the cosine 
law as given in Eq. (6). Hence, substitution of 
Eq. (6) into Eq. (4) gives the constant A in 
terms of y and e„ as 


' * / 

The true situation is that the boundary layer 

mass flow undergoes an inviscid expansion with its 
initial conditions specified by a viscous profile. 
Hence, Eq. (7) is used in Eq. (9) to give 



A = 



2 


TT 



( 8 ) 


where the approximation sin e ~(2e 00 M 
s i n ( tt / 2 e/e*,) is used to simplify the evaluation 
of the integral. This approximation is valid 
unless the integrand becomes negligibly small. 

Next, the angle e 0 and the parameter 
B are estimated. Following Simons* analysis, 
the boundary layer mass flow may be related to 
the mass flow in the plume that is diverging at 
angles greater than e 0 to yield 


The above equation may be simplified somewhat by 
- 39 ^ 

assuming e + 0 (p > 1) to give a simple 
quadratic equation for the parameter e, 
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It should be noted that Eqs. (10) to (12) reduce 
to that of Simons upon application of his simpli- 
fying assumptions for the case of small boundary 
1 ayer nozzle flow. 

Calculations of Nozzle Flow Conditions 


The effects of the displacement thickness 
<s* in reducing the actual area ratio of the 
nozzle are accounted for in this analysis. 
Reasonable estimates of the values of both the 
displacement thickness and also the boundary layer 
thickness may be obtained through an application 
of the Cohen-Reshotko method (Ref. 7). The corre- 
sponding reduced Mach number at the nozzle exit is 
then calculated from a one-dimensional, isentropic, 
ideal flow analysis based on (A e /A*) rec |. Simi- 
larly, standard ratios in terms of y and M 0 
are used to estimate the gas exit velocity and the 
stagnation pressure. The exit velocity is depend- 
ent upon gas type (y), stagnation temperature (T 0 ), 
and nozzle geometry [(M e ) and (A e /A*) rec j] . 


U = 

e e 


1 - 


- 1 ) 


.1 -1/2 

R 


y MW T o 


+ 1/2 

(14) 


Also, for any given mass flow rate m and stagna- 
tion temperature T 0 , the stagnation pressure is 


P 


o 




(15) 


Lastly, the limiting turning angle of the flow 
e„, is defined as the difference between the maxi- 
mum Prandtl-Meyer angle (M e «) and the P-M angle 
for the actual M e plus the half angle of the 
nozzle. Rather, 


Number column density (NCD) is the number of 
molecules per unit area seen by an observer along 
a specified viewing direction, or 1 i ne-of-sight 
(LOS). Column densities can be determined analyt- 
ically by integrating the plume number density 
along a LOS from the observer location to infinity. 
A realistic operational scenario is illustrated in 
Fig. 4. The observer (from which the LOS origin- 
ates) is located at the center of the upper boom 
of the dual keel reference configuration at a 
point 50 m above the transverse boom. The resis- 
tojet source is located at the end of a truss, 
perpendicular to and 35 m away from, the trans- 
verse boom. The observer and the resistojet are 
in the same spatial plane; consequently, this 
represents a worst case situation. 

The resistojet thruster is assumed to operate 
under typical conditions, e.g., using Station 
waste fluids/gases as propellent. The Space 
Shuttle's fuel cells produce waste water which can 
be scavenged and used on the Station for propul- 
sion. Therefore, column density calculations for 
a resistojet using steam as the propellant under 
both hot (1000 °C) and cold (300 °C) flow condi- 
tions are of interest and will be presented. 

Also, the Environmental Control /Life Support 
System (ECLSS), Material Technology Labs (MTL) , 
and Attached Payloads ( AP ) are sources of waste 
gases which will be collected by a Space Station 
Waste Fluid Management System and are available 
for use by the resistojet propulsion module as a 
mixture of gases. An assessment of the column 
densities produced by a typical mixture of helium, 
hydrogen, nitrogen, and argon will also be given. 

Apparatus 


Instrumentation 


e ~ = ''max " v + * (16) 


In this simplified model, then, the included 
angles 2e 0 and 20 * encompass the inviscid core 
and boundary layer regions of the plume, respec- 
tively. There is no provision for flow going 
beyond the limiting angle e^. 

Calculation of Plume Number Column Densities 

The induced environment produced by a resisto- 
jet propulsion module can be analyzed to predict 
possible contamination which may affect Space 
Station sensitive surfaces and attached payload 
optical experiments. The parameters associated 
with quantifying the induced environment include: 
(1) number column density (NCD) - the molecular 
number density integrated along a given line-of- 
sight; (2) backflow - flow turning through angles 
greater than 90° from the nozzle centerline; (3) 
return flux - contaminant molecules impinging on 
a surface due to intermolecular collisions with 
ambient molecules; and (4) particle emission. In 
order to evaluate the impact of these parameters, 
the exhaust plume density field must first be 
defined. This report presents only a number 
column density assessment for a Space Station 
resistojet. A more complete assessment of the 
impact of resistojets on Station science and tech- 
nology activities can be found in Ref. 8. 


Quartz Crystal Microbalances (QCMs) have been 
used in similar investigations as a means of map- 
ping exhaust flow fields. 9 " 14 Basically, the QCM 
used in this study consists of a matched pair of 
precision 15 MHz quartz crystals. The optically 
polished crystals are cryogenical ly cooled to tem- 
peratures sufficient to collect mass on one of the 
crystals (sensor crystal). The other crystal 
serves as a reference. The QCM translates a 
change in the beat frequency between the two crys- 
tals to mass loading. Sensitivity of the device 
is 1 . 56x10" 9 g/cm 2 Hz. Pertinent details of the 
QCM are given in Table 1 while a close-up of the 
device is shown in Fig. 5. 

The quartz crystals are seated inside a gold- 
plated, copper heat sink as shown in the photo. 
Only the sensor crystal is exposed to mass flow. 
The QCM heat sink is cooled to -196 °C (77 K) with 
a continuous flow of liquid nitrogen (LN 2 ) through 
the aluminum mounting bracket. A thermal break 
made of laminated plastic (G-10, not shown in the 
figure) protects the electronics housing of the 
QCM from the extreme temperature of the heat sink. 
The crystals may be set to a desired temperature 
with a control unit. The available set tempera- 
tures range from near heat sink level (as cold as 
possible) to a temperature 80° above that of the 
heat sink (using the QCM internal heater). The 
output of the QCM is sent to a frequency counter. 
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Resistojet Thruster 


Procedure 


Also shown in Fig. 5 is the laboratory 
resistojet used in this study. The main body of 
the resistojet is made of Inconel* while the 
internal heating element is grain-stabilized plat- 
inum. Although it is possible to supply approxi- 
mately 150 W of power to the heater, the mass 
flux measurement experiments described here-in 
consist solely of cold flow thruster operation. 
Carbon dioxide was chosen as the propellant 
because of its potential usage for Space Station 
propulsion and because the gas can condense at 
the LN? temperatures. The mass flow rate was 
0.20 g/sec at 298 K. This corresponds to a thrust 
level of 0.12 N (28 mlbf). Further information 
on the thruster characteristics and operating 
conditions can be found in the first portion of 
Table 2. 

Figure 5 also points out an obvious diffi- 
culty with the use of intrusive probe techniques 
for local, near field measurements of the resisto- 
jet exhaust. Quite often, the extremely small 
nozzle geometry precludes the use of many probe 
types in the vicinity of the thruster exit. The 
QCM and mounting bracket did cause a disturbance 
to the exhaust flow at some measurement locations; 
however, because the heat sink and bracket sur- 
faces were maintained at the saturation tempera- 
ture of LN?, it can be assumed that the greater 
portion of CO? mass flux hitting the surfaces con- 
densed on them, thereby reducing the disturbance. 

Test Facility and Experimental Set-Up 

The mass flux measurements were conducted in 
the 4.6 m (15 ft) diameter vacuum facility at NASA 
Lewis. The vacuum tank has a 4.2 m inner working 
diameter and is 19 m long. The tank was chosen 
for its size and pumping capacity. It is equipped 
with a pumping train consisting of 20 oil diffu- 
sion pumps, backed by four rotary blowers and four 
rotary piston vacuum pumps. The tank is also 
equipped with a LN?-cooled condensor which lines 
the major portion of the inner tank wall for added 
pumping capability. For the resistojet mass flow 
the specified earlier, the tank maintained a back- 
ground pressure of 2.7xl0“ 3 Pa (2x10” 5 tor r) . 

The cryogenically cooled QCM and the resisto- 
jet were mounted on radial positioning arms as 
shown in Fig. 6. The two arms were rigidly attached 
at the end of push/pull rods that extended the 
assembly into the main tank area through a 0.9 m 
(36 in.) diameter port entrance as shown in the 
figure. A retractable aluminum tray supported the 
two axial positioning rods on a moveable test 
cart. The arms could rotate radially through a 
90° arc (±45° with respect true vertical); simi- 
larly, the push/pull rods could be positioned 
axially to permit a mass flux mapping of one 
quadrant of the assumed axisymmetric exhaust. 
Electronic inclinometers attached at the end of 
the push/pull rods were used to mark the location 
of the arms with respect to the true vertical 
position. Because the tank LN 2 baffles were in 
use during the experiments, the sensitive elec- 
tronics of the QCM and inclinometers were kept 
operational by wrapping them with resistive heat- 
ing tape. 


Prior to each mass flux experiment, the posi- 
tions of the QCM and the thruster were calibrated 
for the region of exhaust to be surveyed. Addi- 
tionally, cooldown of the large vacuum tank 
baffles was initiated at least 15 hr before each 
test. The temperatures of several critical baffle 
locations were monitered throughout the duration 
of the test period to ensure a near-equilibrium 
tank background condition. 

The experimental assembly was installed in 
the 0.9 m diameter bell jar entrance to the tank. 

The QCM sensor was internally heated to +90 °C for 
a period of about 20 min to clean the sensor 
crystal of possible contaminants. Once this was 
completed, the axial positioning rods holding the 
QCM and thruster were carefully extended into the 
main tank area and cooldown of the QCM was begun. 
This cooldown generally took 1 to 1-1/2 hr. The 
QCM heat sink was cooled to -196 °C and then main- 
tained at this temperature with a continuous flow 
of LN 2 through the aluminum mounting bracket. A 
150 liter (40 gal) dewar supplied the LN 2 at 140 
to 310 KPa (20 to 45 psia). The QCM sensor was 
set to a desired temperature with the control unit. 
Once the sensor and the heat sink had stablized at 
desired temperatures, background mass flux data 
(e.g., no thruster flow) were taken for a period 
of 30 min or longer. Typically, the background 
mass flux decreased slightly during the tes*- 
period so it was rechecked periodically. Th£ back- 
ground mass flux was usually 2 to 3 percent of the 
centerline flux during thruster operation. 

Back flux Measurements 

Back flux measurements were conducted to find 
the percentage of mass flow turning through angles 
greater than 90° from the thruster centerline and 
also the limiting turning angle. The mounting 
arrangement for the QCM on its positioning arm did 
not permit in-situ rotation of the QCM sensor into 
the flow streamlines. However, it was possible to 
orient the QCM in three different positions: (1) 

with the plane of the QCM sensor crystal perpendic- 
ular to the thruster centerline; (2) with the sen- 
sor crystal plane at 135° from the thruster center- 
line; and (3) with the sensor plane 180° to the 
thruster exit plane to determine back-scatter 
effects. For the back flux measurements, the QCM 
temperature control was set so that the sensor 
reached its lowest possible temperature 
(-189.6 *C) to ensure condensation of the CO 2 . 
Deposition rates were measured with the QCM 
oriented in all three of the above-mentioned 
positions. 

Forward Flux Measurements 

The QCM used in this study had a limited 
dynamic loading capability. High flux rates or 
rapid changes in the flow pattern would cause the 
sensor to be overloaded. This overload would gen- 
erally be characterized by an exponential rise in 
the frequency output to approximately 8 to 10 kHz 
where the signal disappeared. The limitations of 
the QCM dynamic response and also the physcial 
limitations imposed by the measurement apparatus 
precluded the measurement of CO 2 deposition rates 
in the forward flux region. However, the CO 2 
propellant did contain a small percentage of water 
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(220 parts per million as measured by the dew 
point method). In order to get some estimate of 
the forward flux density distribution, the H 2 O 
was used as a tracer. Inherent to this measure- 
ment technique are the assumptions that: the 

concentration of H 2 O was uniform throughout the 
propellant; HoO (and only H 2 O) condensed on the 
sensor crystal; and the water molecules expanded 
in direct correspondence to the CO 2 , at least in 
the near continuum region of the measurements. 

The QCM sensor temperature was set at -175 °C 
so that it would condense the H 2 O mass flux but 
not the CO 2 . The QCM control unit was able to 
maintain this temperature to within ±0.1 °C for 
most measurement locations. Mass flux surveys 
were taken with the water tracer method at axial 
stations of 36.4 and 58.8 cm downstream of the 
thruster. The QCM sensor crystal was oriented 
parallel to and facing the nozzle exit plane. 

The usual test procedure involved moving the 
QCM and thruster to initial measurement locations 
(usually at a radial distance of at least 50 cm 
from one another), establishing a background mass 
flux rate on the sensor, and then beginning pro- 
pellant flow through the thruster. Data were 
taken at 1 min time intervals, recording the OCM 
beat frequency output, the heat sink and sensor 
temperatures, as well as the temperatures of all 
sensitive electronics. Generally, the QCM data 
were taken for a period of at least 5 to 10 min 
at any one location before moving to a new one. 
Because of the limited dynamic loading capability 
of the QCM, it was easier to move the thruster 
with respect to the QCM rather than vice versa so 
this was done where possible. Nonetheless, it 
usually took several minutes to re-establish a 
stable frequency output whenever the QCM was sub- 
jected to any change in the flow pattern. 

Interpretation of Results 

The QCM beat frequency output, when monitored 
over a period of time, could be directly converted 
to a mass flux measurement (g/cm 2 sec). In order 
to convert the measured mass flux to a CO 2 number 
density, the values of two important parameters 
must be obtained. The first parameter is the gas 
velocity at the location of interest. Measure- 
ments of a simulated rocket exhaust plume-*- 5 
have shown that the gas velocity exponentially 
decreases away from the thruster centerline to a 
value which is one half of the exit velocity at 
the limiting turning angle ©«,. A similiar rela- 
tionship for the CO 2 velocity may be obtained, as 

U ( 0 ) = U e exp (-0.0077 e) (17) 

where IL is the gas exit velocity determined 
by Eq. (14) and the angle 0 is in degrees. 

The second parameter is the Capture Coeffi- 
cient, C. It is defined as the fraction of part- 
icles impinging on a surface which are trapped by 
the surface after the collision. 5 Measurements 
of the Capture Coefficients of gases on cryogenic 
surfaces have been made specifically for cryopump- 
ing applications. 16 " 18 Here, the primary concern 
is with a random flux deposition which is depend- 
ent upon an average gas temperature; hence, the 
capture coefficient was assumed to be independent 
of the rate of incidence on the surface. 1 ' The 
measured values of C for both carbon dioxide 


and water as a function of gas temperature may be 
found in Table 3. In the case of directed flux 
deposition, as in the present case, it seems 
likely that not only will some of the incident 
particles reflect off the surface after the col- 
lision, but also some of the particles already 
trapped by the surface may be desorbed by the 
action of impingement. Specific experiments to 
determine the extent of this effect can be per- 
formed; however, due to the preliminary nature of 
the mass flux data presented here, the measured 
values of C at 300 K found in Table 3 will be 
used in the data reduction. 

Results and Discussion 

This section presents results of the analyt- 
ical and experimental investigation of the resisto- 
jet plume. First, the source flow equations 
developed earlier are applied to a laboratory 
resistojet using unheated CO 2 as the propellant. 
Mass flux data taken with a cryogenical ly cooled 
quartz crystal microbalance are converted to 
number density profiles for comparison. The upper 
portion of Table 2 lists the major characteristics 
and operating conditions of the laboratory resisto- 
jet for both analysis and experiment. This sec- 
tion also presents estimates of the level of 
environmental contamination caused by a Space 
Station resistojet during propulsion operations. 
Table 4 lists the major characteristics and oper- 
ating conditions of the Space Station resistojet 
using H 2 O propellant and also a benign mixture 
as the propellant. Comparisons are made between 
predicted number column density levels and 
proposed quiescent requirements. 

Calculations of the Laboratory Resistojet Exhaust 
How Field 

Calculated nozzle exit conditions and plume 
properties of the laboratory thruster are given in 
the lower portion of Table 2. For the flow condi- 
tions specified, the calculated boundary layer 
comprises over 40 percent of the cross-sectional 
area at the nozzle exit. The calculated exit 
velocity is 634 m/sec. It should be mentioned 
that this value compares well with a measured 
value of 627 m/sec, based on thrust measurements 
of this same thruster under similar flow condi- 
tions. iy The resultant exhaust flow field is 
shown in Fig. 7. Plotted in the figure are lines 
of constant CO 2 number density emanating from 
the resistojet source. The dashed lines in the 
figure denote the o calcul ated angles 0 O and 
See at 28 and 90°, respecti vely. The sharp fall- 
off in number density away from the thruster 
centerline demonstrates the relatively benign 
nature of resistojet operation in that its great- 
est impact will be in a very confined region of 
space directly downstream of the thruster. Mass 
flow calculations based on these density contours 
indicate that the inviscid core region of the 
plume comprises 60 percent of the total mass flow 
while 75 percent of the flow is within a symmet- 
rical cone of half angle 35° and 96 percent of 
the flow is contained within 60°. 

Comparison with Preliminary Mass Flux Measurements 

Figure 8 presents a close-up of the above 
calculated flow field along with experimental mass 
flux data taken with the QCM. The temperature of 
the crystals was set as cold as possible so that 
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CO 2 deposition rates were measured. The symbols 
in the figure mark the location of the QCM and the 
flat portion of each symbol denotes the orientation 
of the sensor crystal with respect to the thruster 
exit. The dashed line in Fig. 8 indicates the 
zone where measurements could be obtained. 

An accurate measurement of the background CO? 
number density in the vacuum tank (e.g., no 
thruster flow) is critical to the back flux meas- 
urements. Because the vacuum facility environment 
reflected its previous experience, it was neces- 
sary to take background data immediately prior to 
each data set. For the data displayed in Fig. 8, 
the measured background mass flux rate was consist- 
ently between 0.8xl0 -11 and l.lxlO" 11 g/cm 2 sec. 
This represents an average CO 2 background number 
density of (6±l)xl0 7 molecules/cm 3 . It is felt 
that the background measurements reflect a true 
condition of the vacuum environment (and not an 
inherent drift in the QCM frequency output) 
because: (1) the vacuum tank condition was in, 

to the extent feasibile, a near-equilibrium con- 
dition; (2) both the QCM sensor and heat sink were 
thermally stable; (3) the QCM beat frequency 
increased at a constant rate over a lonq period of 
time; and (4) changes in the QCM output always 
reflected noticeable changes in the tank environ- 
ment (e.g., change in baffle temperature, etc.) 

The back flux measurements with CO? collec- 
tion yielded results as summarized in Table 5. 

With the plane of the QCM sensor facing perpendi- 
cular to the thruster centerline (Test Number 29 
data), no significant increase in the deposition 
rate above the background occurred when flow was 
introduced through the thruster. For all but one 
of the measurement locations, no portion of the 
sensor crystal area was in direct line of sight 
of the edge of the nozzle lip. Conversely, when 
the QCM sensor was off-set 45° from the thruster 
exit plane (Test Number 34 data), the QCM dynamic 
response was exceeded at locations corresponding 
to angles of 87°, 93°, and 97° as measured from 
the thruster centerline. However, at a position 
corresponding to an angle of 99°, the mass flux 
rate decreased sharply to a measureable level 
which was just slightly above the background. It 
should be noted that, at this measurement location, 
54 percent of the crystal area is still within 
direct line of sight of the edge of the nozzle lip. 
These CO? collection data seem to indicate that the 
flow is travelling along a straight path (i.e., 
free molecular flow). Further, and perhaps more 
importantly, the sharp decrease in density between 
97 and 99 suggests a definite boundary to the 
greatest portion of the mass flux (>99 percent), 
but one which is beyond the calculated limiting 
angle, e ro = 90°. This discrepancy can be accounted 
for, in part, by an actual displacement thickness 
larger than that used in the analysis. Further, 
numerical calculations using Direct Simulation 
Monte Carlo (DSMC) techniques have determined that 
flow in the free molecular region of the exhaust 
is dependent upon the flow angle at the nozzle 
lip and also species molecular weight. A 
detailed numerical treatment of the plume, such 
as a combination of the Method of Characteristics 
and the DSMC methods, is necessary to accurately 
model and assess the back flux. 

An attempt was made to measure CO? deposition 
rates in the forward flux region of the exhaust 
(Test No. 30). The QCM was located at a position 


36.3 cm downstream of the thruster exit at a 
radial distance of 50.0 cm. The sensor crystal 
was parallel to and facing the thruster. Because 
of the extremely high flux rates, the dynamic 
response of the QCM was exceeded at this location. 
Consequently, only an estimate of the minimum 
value of the actual CO 2 number density can be 
made. The value, p _> 2xl0 10 molecules/cm 3 , can 
neither support nor refute the calculated flow 
field distribution in the forward flux region. 

Mass Flux Measurements Using H?0 as a Tracer 

In the forward flux region of the resistojet 
exhaust, mass flux data were taken using H?0 as 
a tracer of the CO? propellant. The measurements 
were made along the radii at two specific axial 
stations downstream of the thruster, at x = 36.4 
and 58.8 cm. The CO? mass flux profiles at the 
two stations are presented in Figs. 9(a) and (b). 
The data are plotted as a function of angle off 
centerline. Comparisons of the experimental data 
with calculated density profiles are shown in 
Figs. 10(a) and (b). The data points are curve- 
fit to the calculated profile at a midrange 
measurement location (0 = 27.5°) in Fig. 10(a). 

The same correction factor is used in Fig. 10(b). 
The QCM data, based on a concentration of water at 
220 parts per million in the CO?, were actually 
lower than the predicted values by a factor of 
3.5. This discrepancy could be the result of a 
systematic error in either the experimental tech- 
nique or the analysis. However, the fact that the 
theory and data are in good agreement in terms of 
an overall shape function of the density profile 
gives some confidence to both the experimental 
and analytical techniques. A complete summary of 
the forward flux QCM data is given in Table 6. 

Number Column Density Assessment of a Space Sta- 
tion Resistojet 

Calculated Number Column Densities as a func- 
tion of 1 ine-of-sight angle for the Space Station 
scenario described previously are presented in 
Figs. 11 and 12. The observer zenith (Fig. 4) is 
located at an angle of 0° while the resistojet 
source is at 145 . Also shown in the figures are 
the proposed number column density requirements 
which must not be exceeded for all lines-of-sight. 
(It should be noted that there are no NCD require- 
ments during nonquiescent periods.) 

The column density requirements for infrared 
active molecules, such as H?0, is lO^ 
molecules/cm 2 . The NCDs for a resistojet using 
H?0 propellant remain relatively low in com- 
parison with the requirement until the observer 
begins to look aft. As shown in Fig. 11, the 
NCDs rise steadily until they exceed the proposed 
limit at about 45 aft of zenith, and then begin 
an exponential rise toward infinity for the LOS 
intersecting the resistojet source. Although not 
shown in the figure, the NCDs will drop rapidly 
once past the source to a relatively low level. 

It should be mentioned that lines of sight greater 
than 90° are presumed to be seldom used viewing 
directions for sensitive optical instruments 
located on the astronomical platform because the 
instruments will "see" a greater NCD due to the 
natural atmosphere than that imposed by the 
thruster. Further, 1 ine-of-sight angles which 
pass near the Space Station structure will tend 
to be avoided by sensitive optical instruments. 
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Nonetheless, the obvious disparity between calcu- 
lated NCO values and the proposed limits indicates 
that a resistojet should not use H 2 O propellant 
during quiescent periods if observations are to be 
made along LOS angles between 45 and 150°. It 
should be mentioned here that the solid angle sub- 
tended by this restricted viewing region represents 
only 20 percent of the total 4ir steradian viewing 
space. 

For the case of a resistojet operating on a 
benign gas mixture (Fig. 12), the NCD curve shows 
the same general shape function with LOS angle as 
the H 2 O resistojet. But, since none of the molec- 
ules in the gas mixture are infrared active, the 
limits for molecules of this type are less strin- 
gent, at 10-*- ^ molecules/cm^ . The figure presents 
NCD curves for three specific operational thrust 
levels. For a thrust level of 0.4 N (100 mlbf), 
the quiescent requirements on NCD are not exceeded 
until the observer views along LOS angles which 
pass near, or look into, the resistojet source. 

Concluding Remarks 

The resistojet thruster is currently base- 
lined for Space Station auxiliary propulsion needs. 
Definition of the resistojet flow field is 
required in order to assess the exhaust plume 
impacts to Station science and technology activi- 
ties. This paper has presented initial efforts 
in an analytical and experimental investigation 
of resistojet plumes. Firstly, a modified version 
of Simons' source flow technique for the predic- 
tion of exhaust flow fields was developed for the 
resistojet thruster. The increased effects of the 
relatively large nozzle boundary layer character- 
istic of resistojet nozzle flow were accounted for 
in the analysis. In a parallel effort, mass flux 
measurements of a laboratory resistojet operating 
on unheated CO 2 propellant were obtained using a 
cryogenically cooled quartz crystal microbalance. 
There is qualitative agreement between theory and 
experiment, at least in terms of the overall den- 
sity shape functions in the forward flux region. 
Further, the QCM data indicated that flow does go 
beyond the calculated limiting turning angle (90* 
for the laboratory resistojet) but the amount 
represents only very small fraction of the total 
thruster throughput. 

This paper has also presented the application 
of Simons' method to a Space Station resistojet 
using H 2 O propellant and also a mixture of gases 
as the propellant. The predicted number density 
profiles were used to calculate molecular number 
column density as a function of 1 ine-of-sight 
angle for an observer located at the astronomical 
observation site on the Space Station. For the 
case of a H 2 O resistojet, calculated NCD values 
rise above the proposed requirements for infrared 
(IR) species when the observer looks between 45 
and 150° aft of zenith. This restricted viewing 
region represents only 20 percent of the total 4ir 
steradian viewing space. However, it is not 
recommended that H^O propellant be used during 
quiescent periods if these viewing angles are 
required by an observer. In the case of a 
resistojet using a benign (e.g., noninfrared 
active) propellant mixture, quiescent limits on 
NCD are not exceeded until the observer looks at 
LOS angles which pass near, or look into, the 
resistojet source. 
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TABLE 1. - CHARACTERISTICS OF THE QUARTZ CRYSTAL MICROBALANCE (QCM) 


Crystal frequency, MHz 

Mass sensitivity, g/cm^ Hz . . . 
Operational temperature range, *C 
Beat frequency output, KHz . . . 
Heat sink dimensions, cm . . . . 
Crystal diameter (exposed), cm . 
Crystal area (exposed), cm 2 . . . 


!*.!!!!!!*.*.! i.56xio- 5 hz 
-100 to -190 (cryogenically cooled) 
1 to 8 


2.92 diam. x 3.94 length 

0.851 

0.569 


TABLE 2. - RESISTOJET NOZZLE GEOMETRY/OPERATING CONDITIONS 
AND CALCULATED NOZZLE FLOW/PLUME PROPERTIES 


Nozzle Geometry 


Throat radius 
Exit radius 

R* = 4.06x10-4 m 
R e = 3.81xlO" 3 m 

Half angle 

* = 20' 

Area ratio 

i 

A e /A* = 88 

Reduced area ratio 3 

(A e /A*)red > « 

Thruster Operating Conditions 


Propellant 

C0 2 

Molecular weight 

MW = 44.01 

Specific heat ratio i 

Y = 1.288 

Viscosity (0 298 K) 

v = 15.2x10”° N sec/nr 

Mass flow rate 

m = 0.20 q/sec 

Stagnation temperature 

T 0 = 298 K 

Calculated Nozzle Exit Conditions 


Stagnation pressure 

P = 1.39xl0 5 N/m 2 
0 

Exit Mach number 

M = 5.24 
e 

Exit velocity of the gas 

U e = 634 m/sec 

Limiting gas velocity 

= 710 m/sec 

Average limiting gas velocity 
Boundary layer thickness 

= 532 m/sec (a = 0.75) 
a = 7.87xlO“ 4 m 

Displacement thickness 

a* = 5.33xl0- 4 m 

Plume Properties 


Normalization constant 

A = 1.41 

Parameter beta 

6 = 4.3 

Angle e 0 

e 0 = 28' 

Limiting angle 

e. = 90“ 


deduction due to displacement thickness. 










TABLE 5. - SUMMARY OF MASS FLUX MEASUREMENTS WITH CO 2 COLLECTION 


Test 

number 

QCM/thruster 

orientation 

Angle off 
centerline, 
degree 

Axial 

distance, 

cm 

Radial 

distance, 

cm 

Percent of 
sensor area 
exposed, 3 
percent 

Measured 
mass flux, 
g/cm 3 sec 

CO 2 number 
density, c,d 
molecules/cm 3 

29 

O 

91* 

-0.3 

13.0 

0 

l.OxlO" 11 

7.4xl0 7 


90* 

-0.3 

50.0 

0 

CM 

1 

O 

X 

co 

5.5xl0 7 


fS 

117° 

-6.6 

13.0 

25 

l.lxlO- 11 

0 

00 

t". 


y 

98° 

-6.6 

50.0 

0 

1 

O 

O 

7.4xl0 7 


Q 






>2xl0 10 

30 


62* 

+36.3 

67.1 

55 

>3.7xlO” 8 









34 


b 87* 

+1.3 

25.1 

60 

>4.6xl0~ 8 

>3xl0 10 


n 

b 93° 

-1.3 

25.1 

56 

>6.4xl0“ 8 

>4xl0 10 


N 

b 93° 

-1.3 

25.1 

56 

>4.3xl0~ 8 

0 

0 

X 

^1 



b 93‘ 

-1.3 

25.1 

56 

>5.0xl0~ 8 

>4xl0 10 



b 97° 

-3.2 

25.1 

55 

>7.8xl0" 8 

O'! 

X 

0 

0 



99* 

-3.8 

25.1 

54 

0 

1 

0 

X 

l.OxlO 8 


a The percentage of total sensor crystal area which is in a direct 1 ine-of-sight of the thruster exit (forward 
flux) or nozzle lip (back flux); e.g., this value gives a measure of the shadowing effect imposed by the 
physical surfaces of the QCM heat sink. 

“QCM overloaded at this location; only an estimate of the absolute minimum CO 2 number density can be made. 
“The background CO 2 number density for these tests was (6*l)xl0' molecules/crrP. 

“Tank pressure was 2. 7x10“ 3 Pa. 


TABLE 6. - SUMMARY OF MASS FLUX MEASUREMENTS WITH H?0 COLLECTION 


[Sensor temperature was -175.1 *C.] 


Test 

number 

Axial 

distance, 

cm 

Angle off 
centerline, 
deg 

Radial 

distance, 

cm 

Percent of 
sensor area 
exposed, 3 
percent 

h 2 o 

measured 
mass flux, 
q/sec crrr 

Velocity 
U(e) P 
cm/s ec 

H£0 

number 
densi ty, c 
number/cm 3 

C0 2 

number 
density, ~ 
number/cnr 

32 

36.3 

61.3 

67.1 

50 

9.9xlO~ U 

39 500 

9.3xl0 7 

4.2x10** 



61.3 

67.1 

50 

7.8xl0 -11 

39 500 

7.3xl0 7 

3.3x10** 



54.0 

49.8 

53 

9.4xl0 -11 

41 800 

8.3xl0 7 

3.8x10** 



46.0 

38.4 

59 

5.7x10”** 

44 500 

4.8xl0 7 

2.2x10** 



40.7 

31.0 

64 

1.6X10- 10 

46 300 

1.3xl0 8 

5.9x10** 



35.8 

26.4 

68 

2.7x10”*° 

48 100 

2.1xl0 8 

9.5x10** 



31,8 

22.4 

72 

6. 1x10“ 10 

49 600 

4.6xl0 8 

2.1xl0 12 



27.5 

19.3 

75 

l.lxlO" 9 

51 300 

7.7xl0 8 

3.5xl0 12 



22.6 

15.2 

80 

1.5xl0“ 9 

53 300 

l.OxlO 9 

4.5xl0 12 



16,7 

10.9 

85 

2.4xl0~ 9 

55 700 

1.6xl0 9 

7.3xl0* 2 



16.7 

10.9 

85 

2.8xl0~ 9 

55 700 

1.8xl0 9 

8.2xl0 12 



5.2 

3.3 

96 

3.3xl0” 9 

60 900 

2.0xl0 9 

9.1xl0 12 

31 

58.8 

37.6 

45.2 

66 

1.7x10“*° 

47 500 

1.3xl0 8 

5.9x10** 



34.6 

40.6 

69 

0 

1 

0 

LT> 

CM 

48 600 

1.9xl0 8 

8.6x10** 



26.6 

29.5 

76 

5.4x10”*° 

51 700 

3.9xl0 8 

1.8xl0 12 



22.7 

24.6 

79 

5.9xl0" 10 

53 200 

4.1xl0 8 

1.9xl0 12 



18.9 

20.1 

83 

6.8x10”*° 

54 800 

4.6xl0 8 

2.1xl0 12 


a The percentage of total crystal area which is in a direct 1 ine-of-sight of the thruster exit (forward 
flux) or nozzle lip (back flux); e.g., this value qives a measure of the shadowing effect imposed by the 
physical surfaces of the QCM heat sink. 

“Gas velocity based on Eq. (15) and U e = 634 m/sec. 

“Tank background density (no thruster flow) was 6xl0 7 molecules/cm 3 . 

“Based on measured water concentration of 220 parts per million in the CO 2 . 






NOZZLE 


POWER SOURCE 


HEAT EXCHANGER 


PROPELLANT 


Figure 1. - Schematic of a resistojet thruster. 


Figure 2. - Engineering model resistojet developed for Space 
Station applications. 
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Figure 6. - Experimental apparatus for mass flux measurements. 






Figure 7. - Calculated exhaust flow field of a C0 2 resistojet. 
Propellant mass flow rate, 0.20 g/sec; T q = 298 K. 
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Figure 8. - Comparison of experimental mass flux data (C0 9 
collection) with source flow predictions. z 
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(a) Measured mass flux as a function of angle off 

CENTERLINE (^0 TRACER METHOD). AXIAL STATION 

IS 36. M CM DOWNSTREAM OF THE THRUSTER EXIT 
PLANE. 

Figure 9. - Measured H ? 0 mass flux data in the forward 
flux region. 



(b) Measured mass flux as a function of angle off 
centerline (H 2 0 tracer method). Axial station 

is 58.8 CM downstream of the thruster exit 

PLANE. 
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Figure 9. - Concluded. 
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(a) Calculated versus measured C0 2 number density as 

a function of angle off centerline. Axial station 
is 36.4 CM downstream of the thruster exit plane. 
Data are curve-fit at 0 = 27.5°. 


Figure 10. - Comparison of experimental data with analysis. 
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(b) Calculated versus measured Co 2 number density as a 

function of angle off centerline. Axial station 
is 58.8 cm downstream of the thruster exit plane. 
Data are curve-fit using the same correction factor 
as the previous plot. 


Figure 10. - Concluded. 
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Figure 11. - Calculated number column density as a 

FUNCTION OF LINE-OF-SIGHT ANGLE FOR A RESISTOJET 
THRUSTER USING ^0 PROPELLANT. THE OBSERVER 

ZENITH IS AT A l_-Q-$ ANGLE OF 0° WHILE THE RESIS- 
TOJET SOURCE IS AT 145°. The RESISTOJET THRUST 
LEVEL IS 220 MN (50 MLBp) . 



Figure 12. - Calculated number column density as a function 

OF LINE-OF-SIGHT ANGLE FOR A RESISTOJET THRUSTER OPERATING 
ON A BENIGN PROPELLANT MIXTURE OF W 2 ' N 2' HE ' AND AR ’ THE 
OBSERVER ZENITH IS AT A L“0-$ ANGLE OF 0° WHILE THE RESIS- 
TOJET SOURCE IS AT 145°. 
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